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O , Hyperspherical close-coupling method is used to calculate the elastic, positronium 

0\ ; 

O \ formation and excitation cross sections for positron collisions with atomic hydrogen at 

a^ : 

each inelastic threshold are also analyzed. The adiabatic hyperspherical potential curves 

Q , are used to identify the nature of these resonances. 
> ■ 

^ ■ 1. Introduction 

KA , 

■ The scattering between a positron and an atomic hydrogen is one of the simplest- 

looking quantum mechanical three-body problems. Yet it poses major challenge for 
theorists as well as for experimentalists over the years. The difficulty of formulating a 
satisfactory theoretical approach which accounts for all the possible inelastic transitions, 
including excitation, positronium formation, and ionization processes in a broad energy 
range is responsible for the slow progress in this field. 

In the past, positron- hydrogen atom scattering has been treated theoretically using 
the standard methods generalized from electron-atom scattering, such as the close cou- 



pling method [1-6], the R-matrix method [7], and the variational method [8-9]. There 
are some successes, especially in the lower energy region where few channels are open. 
As the collision energy increases the number of open inelastic channels increases rapidly 
and there are few results in the higher energy region where all the possible final channels 
are addressed. 

An alternative approach, developed in the last few years for treating positron- 
hydrogen atom scattering, is the hyperspherical close coupling (HSCC) method. The 
HSCC has been initially developed for two-electron atomic systems [10], but has since 
been generalized to any three-body systems recently by Igarashi and Toshima [11], and 
by Zhou and Lin [12]. The basic formulation is similar to the earlier work of Archer et al 
[13] but the recent success of the HSCC method owes much to the improvement in the 
numerical accuracy which makes the detailed studies of this collision system possible. 

The HSCC method has been used to study in details the J=0 and J=l partial waves 
in e^-l- H collisions in the energy region below the H(n=4) threshold [14]. Partial cross 
sections to each final state and the resonances below and near each inelastic threshold 
were examined. There are few inelastic cross sections available for comparison in the 
higher energy region except some recent ones from the close coupling method [3] . For the 
resonance parameters, calculations based on the complex coordinate rotation method 
[15,16] have been carried out for a number of states and they compare well with the 
results from the HSCC method [14] . The HSCC has also been used to show that there 
are no resonances in the positronium formation channel at energies above the ionization 
threshold [17]. This is a subject of great controversies [4,6,11,18-21]. The HSCC method 
also has been applied to calculate and compare resonances in the collision systems: 
e~-|-H, e^-l-H and e~-|-Ps [22]. In carrying out these calculations, the same set of codes 
are used since the HSCC method has been developed for the general three-body systems. 
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In this paper we report new results for the 3=2 partial wave for the positron- 
hydrogen atom scattering using the HSCC method. The HSCC method is briefly re- 
viewed in Section 2 and recent numerical implementations are addressed. The results 
for the calculated inelastic scattering cross sections and resonance parameters are pre- 
sented in Section 3. Adiabatic hyperspherical potential curves are used to assist the 
interpretation of the calculated resonances. A final summary is given in Section 4. 

2. The Hyperspherical Close Coupling Method 

First we define the Jacobi coordinates of the three particles. There are three pos- 
sible choices of the Jacobi coordinates, as depicted in Fig. 1. From each set of Jacobi 
coordinates, pi and p2, a set of hyperspherical coordinates can be defined. In hyper- 
spherical coordinates, the wavefunctions are shown to satisfy 

(-^ + Had - 2,iE)^{R, = 0, (1) 

where R is the hyperradius, is the hyperangle and Vt denotes collectively the four 
orientation angles of vectors pi and p2- The adiabatic Hamiltonian Had is the total 
Hamiltonian of the system in the center-of-mass frame evaluated at a constant value of 
the hyperradius and is a scaled mass. 

In the HSCC approach, the hyperradius is divided into two regions at R=i?o- In 
the inner region R< Rq where all the three particles interact with each other, the 
wavefunctions are to be expressed in terms of hyperspherical coordinates. In the outer 
region, the three-body system breaks into a single particle and a pair of particles in 
some bound states. In this outer region, the wave functions are appropriately expressed 
in terms of Jacobi coordinates for each arrangement. 

In the so-called diabatic-by-sector method, the inner region is further divided into 
many small sectors. Within each sector, the wave function is expanded as 
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where we have expressed the wave functions in the body-frame of the system, with the 
z'-axis chosen to be along the hne connecting from the positron to the proton. The 
y'-axis is chosen to be perpendicular to the plane of the three particles and that the 
electron is always on the +x' side. In (2), the normalized D-function has good-parity 
[23], with the u)^s being the Euler angles. The angle 9 is between the two vectors, pi 
and p2, is the hyperangle and v is the channel index. Ra is chosen at the midpoint of 
the sector, J is the total angular momentum, I is the absolute value of the projection 
of J along the body frame's z' axis and running from to J for (— 1)'^P = 1 states and 
from 1 to J for (— 1)"^P = —1 states, with P being the parity. 

We solve the basis functions ^yi{Ra;B,4>) Ra for a fixed I component on the 
(6', (f)) plane using the finite-element method. The hyperradial function satisfies a set of 
coupled difi'erential equations 

("^ " 4^ " ^^^) ^'-'^^^ + T.^,iMR)FMR) = (3) 

where the coupling matrix elements V's are defined explicitly in Zhou and Lin [12]. 

In the practical implementation of the HSCC method, the set of hyperradial equa- 
tions (3) are integrated within each sector. Starting with the innermost sector, the 
integration is continued until it reaches the boundary of the next sector. At this bound- 
ary, the total wave functions are expanded in terms of basis functions of the next sector 
from which integration within the next sector can be carried out. This procedure is 
continued from small hyperradius to Rq where the wave function is matched to an out- 
side solution expressed in terms of independent electronic coordinates pi,P2- From the 
matching, one can extract the K-matrix and the partial scattering cross sections 
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(j)^ 47r(2J+l) 

where k is the momentum of the incident particle. 

In order to be able to calculate excitation and positronium formation cross sections 
to excited states, the matching radius Rq should be chosen at a relative large value, and 
the hyperspherical basis functions also have to be evaluated at large values of R. In the 
large R region, the basis functions for the lower channels are located near the singular- 
ities of the potential surfaces. We improve the numerical accuracy by introducing two 
new angles which are more appropriate for describing the large R region, as discussed 
in Zhou and Lin [14]. Another complication due to the degenerate hydrogenic excited 
states is that the coupling among the different subshells in a given principal quantum 
number is quite large. To accommodate the interaction among the degenerate hydro- 
genic states, we use the dipole states of Gailitis and Damburg[24] and of Seaton[25]. 
The details of this transformation is also discussed in Ref. [14] . 

3. Results and Discussion 

To obtain results for D wave scattering, the basis functions for the 1=0, 1 and 2 
components are first calculated. The matching radius is chosen at i?o=250.655 a.u. 
and 351 sectors are used in the inner region. Since we are interested in energies below 
the H(n=4) threshold, in the solution of hyperradial equations (3) we include 21 basis 
functions for 1=0, 13 for 1=1 and 9 for 1=2. The basis set can be enlarged to achieve 
higher accuracy but at the expense of more computing time. 

In presenting the results, we first mention that the thresholds, in increasing order, 
for Ps(n=l), H(n=2), Ps(n=2), H(n=3) and H(n=4) are located at 0.5, 0.75, 0.875, 
0.8889 and 0.9375 Ryd, respectively, above the ground state of H. There are resonances 
below each threshold. 
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3.1. Elastic and Ps(n=l) formation cross sections 

In Fig. 2 we show the D-wave elastic and Ps(n=l) formation cross sections for 
positron scattering with atomic hydrogen in the energy region between 0.7 and 0.94 
Ryd. Both cross sections change httle in the whole energy region shown. There are 
resonance structures below each inelastic threshold, but pronounced ones are observed 
only below the H(n=2) threshold. In the figure we also show the close coupling results 
of Mitroy and Ratnavelu[3] where they performed close coupling calculations using six 
atomic states on each center. Their calculation will be denoted as CC(6,6). The results 
from the two calculations are in quite reasonable agreement. 

The agreement in the elastic and Ps(n=l) formation cross sections between the 
two calculations, however, does not imply that the predicted resonance positions are 
identical. Both have calculated one Feshbach resonance below the H(n=2) threshold. 
The position and the width of this resonance are listed in Table 1. Note that the CC(6,6) 
results are quite different from the HSCC results. In similar comparison for the S- and 
P-waves results, the CC(6,6) results also show large discrepancies from the HSCC results 
in the resonance parameters. In those cases, the HSCC results are in good agreement 
with those obtained variationally using the complex coordinate rotation method[15,16]. 
The CC(6,6) results were different. We tend to believe that our results here are more 
accurate than those from the CC(6,6) calculations. 

The HSCC method allows one to analyze qualitatively and semiquantitatively the 
resonance structures. In Fig. 3 we show the three adiabatic potential curves that 
converge to the H(n=2) threshold. Only one curve is attractive which can support 
Feshbach resonances. The position of the calculated resonance is shown as dotted lines 
in the figure. Note that one would expect an infinite series of resonances associated 
with the lowest curve since this potential approaches the threshold as — 0.81/i?^. Such 
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a dipole potential can support an infinite number of resonances, but all the higher ones 
are very close to the threshold and were not examined. 

3.2. Excitation cross sections to the 2s and 2p states 

We show in Fig. 4 the D-wave excitation cross sections to 2s and 2p states. The 
results from the CC(6,6) calculations are also shown. Note that the results from the two 
calculations are again quite close, except in the resonance region. We have examined 
three resonances below the Ps(n=2) threshold. The positions and the widths are listed 
in Table 1 and the results for the first two are compared to the CC(6,6) results. Again 
the discrepancy is quite significant. 

Our calculations show small structures near the H(n=2) threshold. We suspect that 
these structures are due to the relatively small basis functions used in this calculation. 

We also show in Fig. 5 the three adiabatic potential curves that converge to the 
Ps(n=2) threshold. Note that the lowest curve is rather attractive. In fact, it behaves as 
— 17.49/i?^ asymptotically. Note that the second curve in Fig. 5 has a small attractive 
part in the inner region, but the potential well is not strong enough to support a 
resonance. 

3.3. Ps(n=2) formation cross sections 

The cross sections for the positronium formation cross sections to 2s and 2p states 
are shown in Fig. 6. These cross sections are quite small, about 0.5% of the cross 
sections for the elastic and the Ps(n=l) cross sections, and about 10% of the n=2 
excitation cross sections. The resonances are quite pronounced in these small channels. 
The resonance positions and widths below the H(n=3) threshold are studied and the 
results are shown in Table 1. There are no other calculations available for comparison 
in this energy range. 

3.4. Excitation cross sections to n=3 states 
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The excitation cross sections to H(n=3) states have also been calculated. The 
results are shown in Fig. 7. These cross sections are about a factor of two smaller 
than the Ps(n=2) cross sections in the same energy region. The resonances are quite 
pronounced at energies below the H(n=4) threshold and the positions and widths of the 
five lowest resonances are shown in Table 1. 

There are no other calculations available for comparison. The small structures near 
the threshold are likely due to the small number of hyperspherical channels used in the 
present calculation. 

Careful analysis of the resonances at the higher thresholds is more difficult due to 
the numerous avoided crossings. For example, we show in Fig. 8 some of the adiabatic 
curves that converge to the H(n=3) threshold. We also show that the uppermost curve 
which converges to the Ps(n=2) limit has many avoided crossings with the H(n=3) 
curves. 

4. Summary 

In this paper we present D-wave results for the excitation and positronium forma- 
tion cross sections in positron-hydrogen atom collisions at energies below the H(n=4) 
threshold using the hyperspherical close coupling method. There are few calculations 
available for comparison in the higher energy region considered in this paper. Together 
with our previous results for the S- and P-wave, we illustrated that the hyperspher- 
ical close coupling method provides a very powerful method for treating general re- 
arrangement collisions. Calculations including the ionization channels are underway 
which would allow the HSCC method to be extended to the higher energy region. The 
method can be generalized to collisions between positrons and alkali atoms by employing 
a model potential description of the target atom. 
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Table 1. Positions and widths of D- wave resonances below different inelastic thresholds 
for the e^+H collision system. The present results are from the HSCC calculation, while 
the CC(6,6) results are from the close coupling calculation of Mitroy and Ratnavelu [3]. 
The energies and half widths are in units of Rydbergs. 



present CC(6,6) 
E(Ryd) |r(Ryd) E(Ryd) ir(Ryd) 





Below H(N=2) Threshold 




-0.250031 


1.65(-6) -0.25006 


3.9(-4) 




Below Ps(N=2) Threshold 




-0.143429 


1.66(-4) -0.13715 


1.25(-4) 


-0.128529 


6.66(-5) -0.12765 


3.3(-4) 


-0.12503 


6.32(-6) 






Below H(JN=3) Ihreshold 




-0.114529 


4. 79 (-4) 




-0.111630 


7.89(-5) 






Below H(N=4) Threshold 




-0.075582 


6.27(-5) 




-0.069775 


6.76(-5) 




-0.066940 


4.70(-5) 




-0.064103 


1.39(-5) 




-0.063855 


2.74(-5) 





11 



Figure Captions 

Fig. 1. Three sets of Jacobi coordinates for the three particle system. 

Fig. 2. D-wave elastic and Ps(n=l) formation cross sections for positron scattering 
with atomic hydrogen. Elastic cross sections: solid lines, present results; dotted lines, 
from CC(6,6)[3]. Ps(n=l) formation cross sections: long dashed lines, present results; 
dash-dotted lines, CC(6,6). The two results for the Ps(n=l) formation are almost 
indistinguishable. The sharp structures are due to the resonances below each inelastic 
threshold. 

Fig. 3. D-wave adiabatic potential curves that converge to the H(n=2) threshold. 
The position of the lowest Feshbach resonance is indicated by the dotted lines. 

Fig. 4. D-wave excitation cross sections to 2s and 2p states for positron scatter- 
ing with atomic hydrogen. Excitation to 2s: solid line, present result; dotted lines, 
CC(6,6)[3]. Excitation to 2p: long dashed lines, present results; dash-dotted lines, 
CC(6,6). 

Fig. 5. D-wave adiabatic potential curves that converge to the Ps(n=2) threshold. 
The position of the first Feshbach resonance is indicated by the dotted lines. 

Fig. 6. D-wave positronium formation cross sections to 2s and 2p states for positron 
scattering with atomic hydrogen obtained from the present calculation : solid line for 
Ps(2s) and dotted lines for Ps(2p). 

Fig. 7. D-wave excitation cross sections to H(3s), H(3p) and H(3d) states for 
positron scattering with atomic hydrogen. Solid line: excitation to 3s; dotted lines: 
excitation to 3p; dashed lines: excitation to 3d. 

Fig. 8. D-wave adiabatic potential curves that converge to the H(n=3) threshold. 
The uppermost curve converging to Ps(n=2) is also shown. 
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Fig 1, Yan Zhou and C D Lin, Hyperspherical Close-Coupling Calculation of D- wave ... 
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Fig 2, Yan Zhou and C D Lin, Hyperspherical Close-Coupling Calculation of D-wave . 
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Fig 3, Yan Zhou and C D Lin, Hyperspherical Close-Coupling Calculation of D-wave 




Fig 4, Yan Zhou and C D Lin, Hyperspherical Close-Coupling Calculation of D-wave . 




Fig 5, Yan Zhou and C D Lin, Hyperspherical Close-Coupling Calculation of D-wave . 
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Fig 6, Yan Zhou and C D Lin, Hyperspherical Close-Coupling Calculation of D-wave . 
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Fig 7, Yan Zhou and C D Lin, Hyperspherical Close-Coupling Calculation of D-wave . 
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Fig 8, Yan Zhou and C D Lin, Hyperspherical Close-Coupling Calculation of D-wave ... 



